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Spectral Analysis of a Discrete Metastable System 
Driven by Levy Flights 

Toralf Burghoff* and Ilya Pavlyukevichl 


Abstract 

In this paper we consider a finite state time discrete Markov chain that mimics the behaviour of solutions 
of the stochastic differential equation 

Xf(x) = a; - f C/'(Af) ds + sLt, 

Jo 

where [/ is a multi-well potential with n ^ 2 local minima and L = {Lt)t^o is a symmetric a-stable 
Levy process (Levy flights process). We investigate the spectrum of the generator of this Markov chain 
in the limit e ^ 0 and localize the top n eigenvalues Af,..., A^. These eigenvalues turn out to be of the 
same algebraic order C>(£“) and are well separated from the rest of the spectrum by a spectral gap. We 
also determine the limits lim£_>o £~“Af, 1 ^ i ^ n, and show that the corresponding eigenvectors are 
approximately constant over the domains which correspond to the potential wells of U. 

Keywords: metastability; Levy flights; a-stable Levy process; eigenvalues; spectral gap; Markov chain; 
transition times; fractional Laplacian; semiclassical limit. 

AMS Subject Classification: 60J10*, 60J75, 47A75, 47G20, 15B51, 15A18 

1 Introduction. Spectral properties of Gaussian small noise 
diffusions 

Let U e C^(R,R) be a one-dimensional multi-well potential (see, e.g. Fig. [TJ such that \U'{x)\ > 

as a; —> +cx) for some c > 0. Assume that U has exactly n local minima trii, 1 ^ ^ n, as well as n — 1 

local maxima 1, enumerated in increasing order 


—GO = So < mi < Si < m 2 < ■ • ■ < m„ < Sn = GO. (1.1) 

Moreover, we suppose that all these extrema are non-degenerate, i.e. 

U" {mi) >0, 1 ^ ^ n, and U"{Si) <0, — 1, (1.2) 

and denote 

M ■■= {mi,... ,m„}. (1.3) 

The points trii, i = are exponentially stable attractors of a deterministic dynamical system 

generated by the ordinary differential equation xt = —U'{xt)- For any initial position x e (si-i,Si), the 
solution xt{x) does not leave the domain of attraction (si-i,Si) and xt{x) ^ mi as f —> go. 

Although deterministic gradient systems with multiple attractors appear naturally in various applica¬ 
tion problems (e.g. climate modelling mnziEH], physics [32] 1 robotics m , economics m, their evident 
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drawback often consists in an impossibility of transitions between the domains of attraction. To over¬ 
come this limitation, a small noise is often added to the model. Intuitively the perturbed random system 
stays in the vicinity of one of the attractors m; most of the time, making sporadic transitions to the 
other domains of attraction. This type of behaviour is often referred to as metastability. The essential 
characteristics of the small noise system as a metastable hopping process are the life times in the domains 
of attraction, transition probabilities between the wells, the most probable state etc. Naturally, these 
characteristics depend heavily on the type of the noisy perturbation. The case of Brownian perturbations 
is the most well studied. 

For a standard Brownian motion W = {Wt)t^o, consider the SDE 

Xl{x)=x-\ U'{Xl)ds + eWt. (1.4) 

Jo 

Under the conditions on U formulated above, for any x e R, (lEl has a unique strong solution X^ = 
)t&o (see e.g. |14p. Since the classical results by Kramers [32], it is known that the mean life times 
of X® in the potential wells are exponentially long with respect to the parameter e, more precisely of the 
order ^ Considering a generic perturbed 

gradient system on exponentially long time scales of the order r(e) = , /i > 0, one can recover 

different behaviours. In particular on the slowly growing (short) scales, when e''^® is less than the life 
times in all the wells (/r < 2([/(Si±i) — U (mi)), i = 1,... n), the trajectory (X|.(g).()tg[o,i] does not succeed 
in leaving any of the domains of attraction. On longer time scales, the trajectory may leave some of the 
domains of attraction thus making transitions and eventually forming a sublimit distribution. On the 
longest scale, X® converges to its stationary law which is concentrated in the vicinity of the potential’s 
global minimum. For a generic potential, there are n — 1 such critical values which are different and 
can be ordered in decreasing order /r 2 > ps > • ■ • > Mn. These values determine the time scales 
T{e) = which distinguish the sublimit distributions. The values fj,i, i = 2,... ,n, can be calculated 

probabilistically in terms of the heights of the potential barriers of U with the help of the large deviations 
theory by Freidlin and Wentzell (see |21l §VI.6] or a recent exposition |13|1. 

Recall that X® is a strong Markov processes with the infinitesimal generator 

Cef{x) = ^f"ix)-U'ix)fix), /eCo”(R,R). (1.5) 

The existence of multiple sublimit distributions is reflected in the spectral properties of the opera¬ 
tor Le- The generator jCe is a negative definite essentially self-adjoint operator in the weighted space 
Z/^(R, da;). It has a discrete non-positive spectrum 0 ^ —A® ^ —Af > • • • and a corresponding 

sequence of eigenfunctions {4?f which form an orthonormal basis (see |24[ Chapter 3] for the detailed 
exposition). Whereas its top eigenvalue is identically zero, Ag = 0, with the corresponding eigenfunction 
being a constant, it turns out that the next n — 1 eigenvalues are of the order e ''"® where the rates fii 
are obtained above. The principle non-zero eigenvalue \i corresponding to the longest time scale e''^'® 
determines the speed of convergence of the law of Xf to the stationary distribution as t —► oo; further 
eigenvalues recover the speeds of convergence to sublimit distributions. With the help of the Fourier 
method of separation of variables, the functional u^{t,x) = Ea,/(Xf) which is the solution of the Cauchy 
problem 8tU = CeU, Ue(0, x) = f{x), can be represented as 

f dy „ ^ 

ue(t,x) = -^ -;-+ S( <^y) + RS,x). (i.e) 

Jr 

The remainder term can be shown to be negligible compared to the leading n terms of the expansion 
due to the spectral gap property, namely, that there is a constant M not depending on e such that 
A^+i > M for all e e (0,1]. Moreover, the eigenfunctions 4?f, 2 ^ i ^ n, are almost constant over the 
potential wells. 

A rich literature is devoted to the analysis of metastability of the system of the type (11.41) on different 
levels of rigour. A probabilistic characterization of the principle eigenvalue of Ce in a bounded domain 
with the Dirichlet boundary conditions was obtained by Khasminski [28] (see also Friedman |22l Lemma 
1.1]). Schuss et al. |35[l42lH5| determined the shape of the eigenfunction $1 and the asymptotics of the 
eigenvalue A|, especially its subexponential prefactor in terms of the values U"{Si) and U"{mi) with the 
help of formal asymptotic expansions of solutions of second order ordinary differential equations. Buslov, 
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Makarov and Kolokoltsov | 10lllll[30ll31[l34 | established a connection between the top spectrum of the 
diffusion’s generator and the spectrum of the matrix of the inverse mean life times of X® in the potential 
wells, found very exact approximations for the corresponding eigenfunctions and proved the existence 
of an e-independent spectral gap. Bovier, Eckhoff, Gayrard and Klein in [6ll8ll2n| developed a potential 
theoretic approach to metastability of Markov chains and gradient diffusions and obtained very exact 
asymptotics of the top eigenvalues. Metastable behaviour of diffusions in a double-well potential was 
studied by probabilistic methods in [2^291391 . Applications of the spectral theory to simulated annealing 
can be found in to stochastic resonance in |24| , and to analysis of molecular dynamics [sail]. 

Berglund and Gentz [4] studied the behaviour of small noise diffusions in potentials with non-quadratic 
extrema. We refer the reader to a recent review [3] by Berglund for further references on the subject. 

In recent time, equations driven by non-Gaussian noise, especially a-stable Levy processes (Levy 
flights), are being adopted for description and modelling of various real world phenomena (see e.g. [36]). 
Let L = {Lt)t^o be a symmetric a-stable Levy process with the characteristic function 

Ee“-^* _g-c(a)thr^ as (0,2), (1.7) 

where c{a) = aj” dy e (0, oo). Such a process has heavy tails and for convenience the constant 

c{a) is chosen to guarantee the following tail asymptotics: 

lim u“P(|Li| ^ u) = 1. (1.8) 

14^00 


For e > 0, the SDE 

X!{x) = x- ( U'{X!) ds -t eLt (1.9) 

Jo 

possesses a unique strong solution X® = (Xf)t^o which is a strong Markov processes (see [Ill27p. 

In the limit e —> 0, the process X® enjoys the metastable behaviour in the following sense. 

Theorem 1.1 (Theorem 1.1, |27p. Let k = 1,... ,n and x e (sfe_i,Sfe_i). Then the process {X^^^a)t>o, 
Xq = X, converges, as e —> 0, in the sense of finite dimensional distributions to a Markov chain Y = 
{Yt)t>Q, Yq = rrife, on the state space Ai with the stable conservative generator Q = (qij)"j=i given by 


rl\ 1 1 I 

2 1 - mi|“ |sj-m,:|“r 

1 / 1 1 


( 1 . 10 ) 


In other words, there is a time scale T{e) = e “on which the process X® reminds of a finite state 
Markov chain Y. It was also shown in [27], that on slower time scales e p < a, the process X® does 
not leave the potential well where it has started. On faster time scales p > a, one cannot obtain a 
meaningful limit of X® since its life times in the potential wells converge to zero and the process X®yj,a 
persistently jumps between different wells. 

The generator of X® is the integro-differential operator 


T^efix) 


fix + z)-2nx) + f{x-z) 
-|e“(-A)“/V(a;) - U'{x)f{x), f e Co”(R,R), 


( 1 . 11 ) 


and the metastability result of Theorem 11.11 suggests that the top spectrum of TLe should remind of the 
spectrum of the Markov chain Y. One can expect that the top n eigenvalues A® = 0, AI,..., A^ of 
should be closely connected with the eigenvalues of the matrix Q and well separated from the rest of the 
spectrum by a spectral gap. Since the time scale in the metastability result is unique, the eigenvalues 
should be all of the same order 0(8°“). 

Unfortunately not much can be said about the spectrum of X>e. No weighted space is known where 
X>e is self-adjoint and no general spectral theory can be applied in this case. Although it is known that 
the invariant distribution exists and is unique | 33[I41| . explicit formulae for its density can be obtained 
only in a very few particular cases (mainly for the Gauchy process a = 1 and polynomial potentials, 

see [mnsirr^ i 

In this paper we make the first step towards a better understanding of the spectral properties of X® 
by reduction of a jump-diffusion to a finite state discrete time Markov chain and analysing its spectral 
properties in the limit e —> 0. 
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2 A discrete Levy driven system and the main result 


We construct a discrete-time Markov chain on a finite state space that mimic the metastable behaviour 
of the solution V® of (HH) in the limit of small e. The construction is based on the standard Euler 
scheme for SDEs. 

For any 7 > 0 there are a range parameter i? > 0, a time step /i > 0 as well as the spacial mesh 
parameter 5 > 0 such that + h + 5 ^ 'y and such that the following holds true. 

All the local minima of U belong to [—i?, R), —R < mi < m„ < R. With R fixed, let us redefine the 
points So := —R and Sn '■= R for convenience. 

There exists a finite set of points S and a partition Ix = [so,Sn) = consisting of the intervals 

Ix = [<ix,bx) such that max^jg^ \bx — ax\ ^ S/2, each interval R contains only one point from S and 
X e {ax:,bx), M c S, {so,.. . ,Sn} c {ax,bx}xes, and for any x e S\M 


X - hU'ix) i Ix. 


( 2 . 1 ) 


Example 2.1. The set S and the intervals R can be constructed as follows. Let R > ^ v |mi| v |m„[. 
Let 5 < ^ A i minij' |si — mj|. 

First, we demand that all mi G <S and set R. = [mi — |,mi -f |). Choose 0 < h ^ 7/3 such that for 
alH = 1,...,n, all a; e [si_i + |,mi — |] 

X - U'{x)h e [si^i,mi] (2.2) 


and for a; e [mi + |,Si — |] 

X — U'{x)h e [mi,Si]. (2.3) 

For definiteness, let us construct the partition of the interval [mi -f |,Si). Let zi be the solution of 
the equation z — U'{z)h = mi -I- |, 21 > mi -I- |. Decompose the interval 


[mi 


-,Z1 ASi--) 


(2.4) 


into a finite disjoint union of the intervals [axjbj,), — ax < f, denote by x the middle point of each 
such an interval, and add x to S. Then due to the fact that U'{x) > 0 on [mi + |,Si — |], we get 
X — U'{x)h G /nil ■ Denote 22 the solution of the equation 2 — U'{z)h = 21 , 22 > 21 and again decompose 
the interval [ 21,22 a Si — |) a finite disjoint union of the intervals of the maximal length |. Continue 
this procedure and assign the last interval to be (ux , Si] . ■ 

For any a: G <S let y*{x) g <S be the unique element such that 


X — hU'{x) e lySf(x). (2.5) 

Let us denote hj T-. S ^ S the mapping corresponding to the operation i.e. Tx ■= y*{x). The 
condition (HU implies that Tx = y*{x) ^ x for x G S\M whereas Tx = y*{x) = x for x e M. 

On the set S define a discrete time deterministic motion = {Z'^)k^o such that 

4 = x, Z'/= y*{Zi_,) = t'^x, k^l. ( 2 . 6 ) 

For any x e S the sequence {Z°(a;)}fe^o is monotone. 

The deterministic motion Z^ mimics the behaviour of the solutions of the deterministic ordinary 
differential equation xt = —U'ixt). It can also be described as a discrete time Markov chain on S with 
the matrix P° = {p% y)x,yi^s of one step transition probabilities given by 


|l, y = y*{x), 
1 0 , otherwise 


Eventually let us construct a discrete-time Markov chain Z^ = {Z/,)k^o on the state space S which 
mimics the behaviour of the jump-diffusion and can be considered as a random perturbation of Z^. 

Denote y := min{i/: y e S} and y := max{y. y e S} and define the matrix of one-step transition 
probabilities P® = {p%,y)x,yes of Z® as 

pl,y ■=P{x - hU'{x) + Li e ly), y^y,y, (2.8) 


as well as 

Px y := P(a: — hU'{x) + eh° Li ^ by) 


(2.9) 
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(2.10) 


Px,y •= IP(® — hU'{Xk) + £h°^ Ll 5: Uy). 

By construction, the Markov chain gets reflected at the barriers —R and R which mimics the very 
fast return from infinity of the trajectory of X® due to the fast increase of the potential at infinity. 

First we have to study the metastable behaviour of Z® as e —> 0. Since Z® is constructed in such a 
way that it resembles the Levy driven jump diffusion X® one could expect that the Z® demonstrates the 
same metastable behaviour. Indeed, the following analogue of the Theorem 11.11 holds true. 

Theorem 2.2. Let k = 1,... ,n and x e S n[sk-i,Sk]- Then the process (Z f t .)«>n. Zn = x, converges, 

' he^ * 

as e —> 0, in the sense of finite dimensional distributions to a Markov chain Y — (yt)t>o, To = rrifc, on 
the state space M with the generator Q defined in (II.IOII . 

Proof. The proof of this result essentially follows the arguments presented in m for the case of a one 
dimensional jump diffusion X®. The arguments are even easier since by construction, S does not contain 
the local maxima of U Furthermore thanks to the reflection condition at ±i? and the finiteness of <S, we 
do not have to consider returns of Z® from infinity. Note that this convergence result does not depend on 
the size or the particular choice of the state space S. For a detailed proof we refer the reader to Chapter 
3 in [^. □ 

The Theorem 12.21 suggests that top spectrum of the generator of the Markov chain Z® should be 
closely related to the spectrum of the matrix Q. Let N — |5| denote the cardinality of the set S. To 
compare the generators of Z® and Y we note that 

P®=e^""’", ( 2 . 11 ) 

so that the matrix P® — Ijv can be viewed as the discrete analogue of an infinitesimal generator. Taking 
into account the time scaling t i—> which appears in the metastability result (Theorem 12.21) . we 

introduce the N x N matrix 

Q":=;^(P"-Iiv) ( 2 . 12 ) 

and denote by o'(Q®) := {A®,..., A^} and o'(Q) := (A)^, ..., A)?} the spectra of Q® and Q, respectively. 
Then the following main theorem holds. 

Theorem 2.3. The spectrum cr(Q®) can he divided into two disjoint parts (Ji(Q®) and cr 2 (Q®) for which 
the following assertions hold: 

(i) o'i(Q®) contains precisely n eigenvalues A®,..., A® such that A® = A^ = 0 and 

limA® = A9 2^i^n; (2.13) 

e^O 

(a) For any seguence {re}e>o such that re —> +QO, ree"+i —> 0 there is £o >0 such that for all e e (0, eo] 
and all i = n + . ,N 


lAfl^r,. (2.14) 

The proof the Theorem [2]3] does not take unto account any a priori information about the eigenvalues 
of Q other than the fact that one of them is 0. However, we strongly believe that in the generic case, e.g. 
when the potential U does not have any intrinsic symmetries, all eigenvalues of Q are real and simple. 
This conjecture is supported by numerous computer simulations of spectra of a matrix Q, when the 
entries of the form (II.IUII are randomly generated. 

Denote by = {5^)xeS, 1 ^ i ^ n, the indicator function (vector) of the set Sj S n [si_i,Si], i.e. 
(5J = 1 for X e Sj and 0 elsewhere. 

Theorem 2.4. Assume that the eigenvalues of Q are real and simple. Let 1 ^ i ^ n and let = 
(^/)®’*)a ,£5 be the right eigenvector of Cf associated with A® and normalized such that maxi^^jcn IV’m’jl = 
1. Let = ('0^’*)j=i be the right eigenvector of Q associated with and normalized such that 

maxi^^jjn = 1- Then the following limit holds: 

71 

mj^l'i/'S’* - ^^ 0, £ ^ 0. (2.15) 

In particular, for i = 1, Af = 0 and = 1. 
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Figure 1: A potential U with the minima rrii = —4.015, m 2 = 0.468, m 3 = 3.966 and the maximasi = —1.034, 
S 2 = 1.921; the range i? = 5. 



lyi, 


1 ■ 10 “ 
5 ■ 10 ® 


-5 ■ 109 


-1 ■ 10 ^ 


Figure 2: The spectrum of the matrix Q® on the complex plane. The top three eigenvalues are located in 
the neighbourhood of the origin and well separated from the rest of the spectrum. 


A V’! 


0.5 


-5 -4 -3 -2 -1 


-0.5 


e,2 


A i’x 


0.5 


1 2 3 4 5 X -5 -4 -3 -2 -1 


-0.5 


123453; 


Figure 3: The eigenvectors and of the matrix are almost constant over the domains of attraction 
of [/. 


Example 2.5. Consider a three-well potential U such that, for simplicity of simulations, U' = +1 outside 
of small neighbourhoods of the extrema, see Fig. [T] and fix a = 1.8. The matrix Q calculated according 
to (II.lot has the eigenvalues = 0, A^ = —0.124, A^ = —0.579, and the corresponding eigenvectors 
- 0*54 ^ 0 * 5.2 _ (—0.629, 0.280,1), 0^^’® = (—0.088,1, —0.245). We construct the state space S 

consisting of A = 203 points from the interval [—5, 5) and a Markov chain with the matrix of transition 
probabilities defined according to (12. 8 t . (12.9t . (12.lot with h = 10/A and e = 10 ®. The spectrum of the 
matrix consists of 203 eigenvalues, see Fig. [2] and the top eigenvalues are Af, A 2 

and A 3 , coincide with the spectrum of Q with a four decimal places precision. The eigenvector 0®’^ is 
constant on S whereas the eigenvectors 0®’^ and 0®’® converge to the values 0^’^, 0^’^. k = 1,2, 3, on 
the parts of S corresponding to the domains of attraction of U, see Fig. [3] 
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3 Eigenvalues. Proof of the Theorem 12.3 

Let Li by a symmetric a-stable random variable with the characteristic function (EZl). It is well-known 
(see [3) that Li has a density Pa{x) on R which can be expanded as a; —> -l-QO as: 

In particular, for any fixed 0 < a < b and h > 0 there is C > 0 such that in the limit e —> 0 

|p(a ^ < 6) - ^ (3.2) 

The next Lemma follows directly form the formulae (12.811 . (12.911 and (12.1011 . 

Lemma 3.1. Let a e (0, 2) and the set S and the partition [—i?, R) = fixed. Then there 

is So > 0 small enough and a constant C = C{h, S, R,a) > 0 such that for all 0 < e < eo the following 
estimates hold: 

(i) for any x e S and any y e S such that y # y*{x) we have 

\pl,y - ^ Ce^“ (3.3) 


where 


2 I Iflj, — a; -I- hU'{x)\°‘ |foj, — a; -I- hU'{x)\'^ 


, y t {v,y,y*{x)}, 


^ 2\hy-x + hU'{x)\‘^' 


^ 2|aj; — a; -I- hU'{x)\°‘ ’ 


y = y, 
y = y-, 


(3.4) 


(ii) for x e S and y = y*{x) we have 


\pI,v*(^) - (1 - f^,y*(^) £°‘h) \ ^ 


where 


fx,y*(x) 2 (’ 


ay*(x) - x + hU'{x)\°‘ - x + hU'{x)\' 


:)■ 


(3.5) 


(3.6) 


3.1 Proof of Theorem 12.31 

In the first step we will analyse the structure of the matrix 

A"(A):=Q"-AIv = ^(P"-Iiv)-AIjv, A e C, (3.7) 

which obviously plays the essential role in the investigation of the spectrum of Q®. 

First we observe that the matrices P® and A®(A) possess a block structure. Indeed, recalling the 
decomposition of the state space <S = <Si u ... u Sn, Si = S n [si_i,Si], i = 1,..., n, as a disjoint union 
of states belonging to different potential wells, we may write 







/A 5 , , 5 , (A) 


All ,5„ \ 

p® = 




, A®(A) = 





\Ps„.Si 


Pll.sJ 


\ A|„_5j 


A|„,s„(A)/ 


The blocks P|^,Sj of P'^ determine the transitions of the Markov chain Z® between the wells Si and Sj. 
Only the blocks A|. 5 .(A), i = 1,... ,n, that include the main diagonal depend on A. 

For the further analysis of the matrix A® (A) it is helpful to study the limit behaviour of its entries 
a® y(A) as e —> 0 . 

The following expansions follow immediately from Lemma l3.1l 
Lemma 3.2. Let A e C, then 

(i) for x,y e S and y f {^x,y*{x)] we have a%^y = dx,y{l- + (!?(e“)); 
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(ii) for X e S\A4 and y = y*{x) we have a® = 0{£ “), more precisely, 

o-l,y*(x)W = where 3 ^,y*( 2 :) =/ 3 ;,;,*( 2;)(1 + ^(e^)); (3.9) 

(Hi) for X G S\Ai we have = 0{£~°‘), more precisely, 

\ + al^,^ where 3®,^ = + 0(e“)); (3.10) 

(iv) for X e A4 we have a% ,^ = 0(1), more precisely, 

a^,x = -3^,x-A where 3®_„, =/:„,,i(l + 0(e“)). (3.11) 

As in Lemma \S. 1\ the bounds hidden in the Landau symbols O are uniform over x,y e S. 

To get a better understanding of the structure of A®(A), let us for example take a detailed look at 
the block ^^(A). If we number the states of <Si in the increasing order by xi,... ,Xn(i), then this 
block has the form 

(3.12) 


^mi,mi A 






mi.":Ar(l) 


_\ -I- n® 


Note that the numbers 3® ,, that appear in the previous lemma are all positive. 

Let us establish the connection between the matrix ^^(P® — In) and the matrix Q of the limit 
Markov chain appearing in the Theorems 1 1.1 1 and 12.21 

Lemma 3.3. Let Q = (q'i,j)"j = i be the generator of the Markov chain Y on the state space M defined 
in (11.101) . Then there exist constants C = C{5, h, R, a) and £o > 0 such that for every 0 < e < eo, every 
I iii,j ^n, i A 3, 




t/eS,- 


and 


y^s. 

Proof. Let for definiteness 1 ^ i < j < n. Then according to (iTSl) . (IT^ and (ITTOI) 

E P-i^y = E ^ 

yeSj 

= +£“3Li g [Sj_i,Sj)) 

= I (7- - -^ - 7- ^-r-) + C’(£“) = <lhj + C>(e“)- 

2V(sj-.i-m,:)“ 


yeSi 


(3.13) 


(3.14) 


(3.15) 


Analogously, for j = n, with the help of (12.1011 . (13.211 and (11.101) we get 

^ E E ph - ^ 

y<^Sn 

^P^mi + e'^hLx g [Sn_i, +co)j 

+ C?(e“) = + 0 (e“), 


y^Sn 


(3.16) 


2 (sj-i — mi)“ 
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(3.17) 


y<^Si 


and for i = j 

E {pk,y - «t) = ^( - 1 + E - hU'imt) + s^hL, e /,)) 

V<^Si 

= + £°hLi ^ [si_i,Si)j 

= -\{-, -^ + 7- ^-T-) + C’(e“) = <ii,i + C(e“)- 


□ 


Now let us come to the investigation of the eigenvalues. We start with a well known result. 
Proposition 3.4. All eigenvalues of Q and Q® have non-positive real parts. 

Proof. This result follows from the Gershgorin circle theorem (see, for example, Theorem 1.11 in |46p. 
Indeed, since qi^i < 0 for all 1 ^ i ^ n, and < 0 for all a; e 5, as well as 

J] Qi.j = -Qi:i and ^ qx^y = -ql^x (3.18) 

j¥=i y^x 

all the Gershgorin circles are located in the negative complex half plane and touch the origin. In particular 
the principal eigenvalues equal to zero, Af = A^ = 0. □ 

The proof of Theorem 12. 3l mainlv consists in the investigation of the asymptotics of the characteristic 
polynomial P^(A) := det(A®(A)) in the limit e —> 0. We will show that, after an appropriate scaling, it 
converges uniformly on a sufficiently large disc to the characteristic polynomial P’^(A) := det(A*^(A)) 
where A*^(A) := Q — AI^. 

Proposition 3.5. There exist constants eo > 0, and C = C{h, 5, a,R) > 1 such that for every 0 < £ < £o 
and A e C 

|(-£“/i)^""'P‘"(A) - P‘^(A)| ^ CdAI + C')"+^£“. (3.19) 

Before we start with the proof of this Proposition let us show how it implies the statements of Theorem 

m 

Proof of Theorem ] 2. 3[ Both statements follow from Proposition 13.51 and the Hurwitz Theorem (see, for 
example. Theorem 1.3.8 in |40p: 

Theorem 3.6 (Hurwitz). Let G c C be a region and let {fk}k^i be a sequence of analytic functions on 
G that converges to a non-zero function f uniformly on every compact subset of G. Then zq e G is a 
zero of f with multiplicity m if and only if there exists a neighbourhood H c G such that, in every ball 
{z e C: |2 — zo| ^ P} c H, each function ft whose index exceeds some bound ko = ko{R) has exactly m 
zeros, counted according to their multiplicities. 

For the proof of (i) fix a closed ball centred around 0 that includes all the eigenvalues of Q. Then 
Proposition |33] ensures the uniform convergence inside of this ball. Moreover, let a sequence re > 0 be 
such that 

lim Tg = +00 and lim£“rg’''^ = 0. (3.20) 

e->0 e^O 

Then, Proposition [33] implies that 

lim sup |(-£“/i)^“"P"(A)-P‘^(A)| = 0 (3.21) 

and that proves the statement (ii). ( Theorem \2.tM ~\ 

The rest of this section is devoted to the proof of Proposition [33] We start of with some preparations. 
Recall that the set S consists of N elements, let 

Sn = S| 5 | := {tt = {'Kx)xcS '■ TT is bijective mapping 5 —> 5} (3.22) 

be the set of all permutations of S, and denote by id the identity permutation. The Leibniz formula for 
the characteristic polynomial P®(A) gives 

P^(A)= ^ sgn(7r)n<..(A). (3.23) 

ttsSat xeS 
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Let us briefly explain the heuristics behind the following analysis of P®. Let the parameters h, 5, R 
and a be fixed. Clearly, in view of the asymptotics of the entries a% y as e —> 0, for any fixed A e C 
we have |P‘^(A)| —► oo as e —> 0. With that in mind, let us omit for a moment the values a® since 
they are of the order 0(1), that is negligible in comparison to e~°‘. Then one can think of P^(A) as of a 
polynomial of the variable e~°‘ with coefficients depending on A and terms of order 0(1). 

The degree of this polynomial does not exceed N — n since only the rows indexed by the minima 
mi,...,mn do not contain the term . Hence we can write 


N — n 

P^(A) ^ ^ cKA)(£“/i)'L (3.24) 

j=o 

By a multiplication by the factor we single out the coefficient c5;_„(A) since all the other 

summands vanish in the limit e —> 0. On the other hand, c%_„{X) contains expressions of the form 
j^Pmi,yy 1 ^ *0 ^ which, as we saw in Lemma 13.31 are approximations of the entries qtj of 
the generator matrix Q. Finally, having that in mind we prove that c^_„(A) approximates P^{X). 

In view of the Leibniz formula (13.2311 we should first turn to the permutations that contribute the 
terms of order to the polynomial P"^. As we see from Lemma 13.21 in') and (iii), such 

permutations are precisely those which satisfy the condition = x ot = y*{x) for x e S\A4. Consider 
the family of permutations 

Ho := {tt = (tvj:)xeS e Sjv : tj'j: = X or tyj: = y*{x) for all x e <S\A4} . (3.25) 

Lemma 3.7. For every i = 1,... ,n, every y e Si there exists a permutation tt e Ho such that 

TVmi = y- (3.26) 

Proof. Fix y e Si. If j/ = mi, then we can choose tt = id. Otherwise there exists a fc = k{y) 5= 1 such that 
((P°)(*^))y^mi = 1, that is the deterministic motion, if started in y, reaches the minimum mi in k steps. 
Recall that the mapping T: S ^ S, T{x) y*{x). Then the equality ((P°)^*^)!/,mi = 1 is equivalent to 

T'^y = mi. (3.27) 

Now we can define a permutation tt : <S ^ <S by 

■= Ty, TTTy = T^y,..., nj.k-iy = P'^y = mi, km, = y (3.28) 

as well as 

=z iizeS\{y,Ty,...,T'^y}, (3.29) 

It is clear that tt e Hq. □ 

Consider the family of vectors 

Tn := {{yi ,..., yn) ■ yi g Si}. (3.30) 

Lemma 3.8. For every {yr ,..., yn) e Tn there exists a permutation tt e Hq such that 

TTm, = yi, 1 ^ ^ n. (3.31) 

Proof. Let TTi, i — 1,... ,n, he the permutations constructed in the proof of Lemma [3.71 Note that each 
of them is non-identical only on the set Si. Hence a composition of TTi, i = 1,..., n, has the desired 
property. □ 

Lemma 3.9. For every (j/i,..., yn) e Tn there exists a permutation tt e Hq such that 

'Xmi=yj, (3.32) 

Proof. Fix {yi,...,yn) e and an associated permutation tt from Lemma |3.8I such that nm, = yi, 
i = 1,..., n. Since no condition on the images of the states mi,..., mn is imposed in the definition of 
Ho, we can permute them among each other to obtain tt e Hq with the required property. □ 
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For the investigation of the determinant of the matrix A® (A) we need to introduce the following 
disjoint decomposition of the set of all permutations Sn- For p = 0,..., N — n let 


(3.33) 


(3.34) 


rip := {tt e Sn '■ there are precisely p states xi,..., Xp e S\A4 

such that t i^q)} for all 1 ^ 5 ^ p}. 

With this definition we can write 

N — n 

S'iv = □ Up. 

p=0 

Also note that the definition of Do in (13.251) coincides with the one given in (13.3311 . For any tt e Sn we 
also define 

5 ^(7r) := {x e iS\A4 : = x}, 

S^(7r) := {x e S\M : tTx = y*{x)}, 

S^{n) := 5\(^5^(7r)|j5^(7r)). 

Clearly these sets are disjoint and S = LJi=i ‘ 5 *( 7 r). 

The following Lemma follows directly from the definition of lip. 

Lemma 3.10. (i) For every tt G Ho we have = M. 

(a) For every 0 ^ p ^ N — n and every tt G lip we have |<S^(7r) u <S^(7r)| = N — n — p. 


(3.35) 


□ 


Let us introduce functions fp: tt 
iS\AI, namely 


that count the number of fixed points of tt G IIp on the set 


/p( 7 r) := \{x G S\M ■■ -Kx = a;}| = |<SA 7 r)|. (3.36) 

By definition of lip, it follows that fp has values in {0,..., A — n — p}. We obtain a trivial equality 

n = n (3.37) 

2165 ^( 77 ) x€5^(7r) 

which allows us to write the characteristic polynomial as 


N — n - 

'(^)= E f E sgn(7r)(-l)*W . Yl (^ + (A-ax..)) 


p=0 TTGlIp 


(3.38) 


^ Y\ 

3 :G*S 3 ( 7 r) 

Using Lemma fii) along with a simple straightforward calculation yields that there are numbers 
A%, 7 r(A) - 0 ( 1 ) such that 

_1 A 1 A' 




(3.39) 


2 :G<S^ (tt) 


a;G«S2(7r) 


As explained before P® can roughly be viewed as a polynomial in the variable {e°‘h)~^ and our aim is to 
single out the coefficient that belongs to the highest power. 

After multiplying the polynomial P® by the factor (e“h)^“" and taking Lemma [3.1UI (il into account, 
we split the sum in (13.3811 up into three parts to get the representation 

n 

(e“h)^""P®(A) = ^ sgn(7r)(-l)f‘>^"> ]^a^.,^,„.(A) 

ttgIIo *=1 

N—n—l n 

+ ^ sgn(7r)(-l)f"(-)( ^ /?f.o,„(A)(£“h)^—')n<,.„,(A) (3.40) 

ttgIIq 1=0 i=l 

N—n N—n—p 

+ ^ 2 sgn(7r)(-l)f-t-)( ^ A%.,(A)(e“h)^-’^-') H 


P = 1 TTGlIp 1 = 0 

Let us continue the investigation of the permutations. 


X^S^{'T 7 ) 
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Lemma 3.11. Fix tt G Ho such that TTm^ e Si for every 1 ^ i ^ n. Then 

sgn(7r)(-l)^“(-) = (-1)^-". (3.41) 

Proof. Fix 1 ^ i ^ n and let TTm^ = yi e Si. From this condition and from the definition of IIo it follows 
that for every x e Si we have e Si. Hence, we can write tt as a product of n disjoint permutation 
cycles 


TT = Cl • ■ • ■ 


(3.42) 


where each cycle C* acts only on Si. Define ki = k(yi) ^ 0 as the number of steps the deterministic 
motion needs to reach if started in yi (see the proof of Lemma 13.71) . Then we can write each cycle C* 
as 




yi ^ Tyi 

THi lUi, 


T^i ly. y.^ 5: 1^ 

ki = 0 . 


(3.43) 


If case ki = 0 then tt acts as identity in Si, i.e. = x for all x e Si. For the consistency of notation in 
this case, however, we do not omit the identity cycle Ci having in mind that the sign of an identity cycle 
is 1. Let |Ci| denote the length of the i-th cycle. Then \^i\ = ki + 1 and 

sgn(7r) = (-l)l«il-i ■ ... • (-l)i«"l-i = (3.44) 

On the other hand, from the definition of fo it follows that 

n n 

/o(7r) = X! (I‘5i| -{ki + l)) = N-n-J]ki (3.45) 

i=l i=l 

which proves the claim. □ 


We constrnct a farther decomposition of the set Hq. Let Sn denote the set of all permntations 
(T: {1 ,..., n} —> {1,... , n} and set 

n((T) := {tt G Ho; TTm- G 5a-j, 1 ^ i ^ n}. (3.46) 

Then it is easy to see that 

Ho = □ n(ff). (3.47) 

Lemma 3.12. Let cr G Then, for every tt g n(cr) 

sgn(^T)(-l)^-" = sgn(7r)(-l)^«'-\ (3.48) 

Proof. Let a = Ci ' • • ■ ' Ci be a representation of cr as a product of 1 ^ L ^ n disjoint cycles. In 
particular, |Ci | + • ■ • + |Ci I = 

From the definition of n(cr) it follows that there exists a collection of states {yi ,..., y„) e In such 
that 


Unti — y^l ^ S(Ti , . . . , — y^'n ^ S(Tn • (3.49) 

Now, every tt G Hq can be written as a composition tt = tTo- o tt where tt G Hq such that TTm; e Si 
and where tTo- only acts on the states yi = TTm^,... ,yn = according to the permutation cr. So let 
(y(yi,... ,yn) denote a cycle representation of tTo-. With the notations introduced in the previous lemmas 
we then can write tt as 

TT = cr(yi, ...,yn) ^ yi ^ ^ ■ ... ■ ^m„ ^ y„ ■ 1 -^ . (3.50) 

Keep in mind that the numbers ki also depend on yi and that one always has T^'yi = iTii. Since 
fo does not take Tni,...,mn into account we have /o( 7 r) = /o( 7 f). Combining that with the equality 
sgn(cr) = sgn(7ro-) and Lemma l3.121 yields 

sgn( 7 r)(-l)^°''^^ = sgn( 7 ra)sgn( 7 r)(-l)^°^’^^ = sgn(cr)(-l)^“’^. (3.51) 

□ 
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Lemma MU allows us to rewrite the first summand in (13.401) as 


n 

E sgn(7r)(-l)'^“("^ ]^a^. 

ttgIIo i = l 

n 

o-^Sn 7ren(o-) i=l 


n 

= Y sgn(o-)(-l)^“" Y 

(T€Sn 7r£n(cr) i=l 


n 

cr^Sn yi^Sfj^ y-n^San * = 1 


n 

^ (_l)V-n ^ sgn(fT)]^ Y O'k.vW- 


a-eSn y€S(j^ 


This representation makes clear why we also included the term (—into the normalization factor 
for P® in Proposition 13.51 

The next lemma provides us with the crucial estimate. 

Lemma 3.13. There exist constants eo > 0 and C = C{h, 5, R, a,n) > 1 such that for every A e C and 
every 0 < e < eo 


I Ti S sgn(o-)]^a^^.(A)| ^ C'((5+|A|)" ^e“. (3.53) 

(TESti i=ly€Scrj^ (T€Sn *=1 

Proof. First we need a generalization of the binomial formula that can be easily derived by a straight¬ 
forward calculation. Fix numbers ci,..., Cn, di,..., d„ e R. Then 

n n 

n(ci + di) = (3.54) 

i=l i=0 

where 

n 

Go=n* Gi= Yi n cfc n ( 3 . 55 ) 

i = l J:6{ii,... ,ii} 

Now, abbreviate the left hand side of (13.5311 with (LHS). We apply (13.5411 with a = Yjyes ~ 

°%rTi ('^) di = a^^_ (A) to obtain 


(LHS)s: ^ 

(TESn 

- s 


i=l \ y€Sa, 


y^Sn 


n-+E E 


n 


Ck 


n 


i=l i=l l€{l,...,n} 




y^Sa 


E n 

(T€Sn 

n—1 

E E E n 

rreSn ^=l 


E “m,,!,(A)-a^^JA) 

yeSa. 


n 




(3.56) 




(*) 

Since each row sum of P® — liv is equal to 0, we have for every x e S 


Y (P" - liv).,. = - 2 E (P" - (3.57) 

yeSi j = lyESj 

3=^i 
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For (T e Sn define F{(t) {1 ^ i ^ n: ai = i} as the set of all the fixed points of cr. Then, 

help of Lemma 13.31 we derive for some constant C > 0 


s n 

cr€Sn 

= s 

(TESn 

= s 

(TESn 




y^Sa 


n 

i=l 

i0F(fT) 

71 

n 


. ^mi,y 




n 

i=l 




j = l ysSj 


ieF(rT) j^i 


i = l 
j¥=i. 


i=l 

i^F(cr) 




^ E 

a-eSn 

cx\n 


y^S(j 
Q!xn-|F(cr)| 


n 


E E ^a^P'^i^y + E 


iEF(cr) j^i 


j = l ysSj 


t = l 


^ n! (nC'£“)" . 

Similarly one can estimate the expression (*) in (13.5611 as follows: 

n IE J■ n 


y^Sy 

k^F{(T) 


/6{l,...,n} 

ltF{<7) 


n I - E E +E ■ n 


fee{ji,...,jj} Syi=ky€S, 
fceF(cr) 


s^k 


ie{l,...,n} 

l<t:F(^) 


^ (nCe^rdQloo + |A|r-L 
We finish the proof by continuing in (13.561) : 

(LHS) ^ n! • (nCe^)" + n! V ( | (nC'£“)'(||Q||oo + |A|)”“‘^ C max ([|Q|oo + |A|)*'£“. 


Assigning C to be the maximum of C and 1 + ||Q||oo yields the result. 

We now complete the proof of Proposition 13.51 
Proof of Provosition VSf5\ First, we combine the equations (13.401) and (13.521) to obtain 

n 

(-£“h)^""P®(A) = Yj sgn('^)n E 

CTEStl *=l2/£‘5cr^ 

iV—n—1 71 

+ (_l)^- E sgn(7r)(-l)^«'-)( 2 /3f,o,.(A)(£“h)^-"-')n W 

Trello ^ = 0 * = 1 

N—71 N — 71 —p 

+ (-!)"'■" E E «gn(7r)(-l)*(-)( 2 A%,.(A)(£“h)"'—*) n <-■ 

p—1 Trellp / —0 xEtS^(7r) 

An application of Lemma 13.131 then yields 

|(-£“/i)^""P®(A) - P‘^(A)| s: C{C + |A|)"-^£“ 

N — 71—1 71 

+ E I E a'^,o,^(a)(£“/i)^— 

Trello 1 = 0 i=l 

N — 71 N — n—p 

+ E E I E n “S.- 

p—1 Trellp l — O xetS^(Tr) 


with the 


(3.58) 


(3.59) 


(3.60) 

□ 


(3.61) 


(3.62) 
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We have to estimate the second and third summand. For this let us use Lemma [3d0l to renumber the 

states of <S^( 7 r) and <S^( 7 r) such that 5 ^( 7 r) = {xi,... and 5 ^( 7 r) = ,a;jv-n-p}. 

Then, from (13.391) it is easy to see that the numbers A%,Tr(''') are polynomials in A of degree | 5 ^( 7 r)| = 
/p( 7 r). Expanding them with respect to A yields for I = 0 

fp(^) N—n—p 

/3§,P,w(A) = n • n (3-63) 

as well as for — n — p 

/ N — n — p \ 

PivA^)= Ya n n i-s-AA]- ( 3 . 64 ) 

lsS/Jl<...<p;sSJV-n-p \ V = 1 u = fp(TT) + l / 

Then, by Lemma 13.21 and Lemma O we can find a constant c = c(a, h, S, R) such that for 0 < e < eo 


fp(^) N — n—p /?(■"■) N N—n—p 

i/3g,p,„(A)i = n n n n 

1^=1 i ^ = / p (7 r ) + l 




U = 1 

/p(7r) 


1 

K,=j^U 


< 


’'=fp(-^)+i 


( n (|A| +iVc))(c)'^-"-^-^'’'’' 


(3.65) 


U=1 


C (|A1 + 
^ (|A| + 


and similarly 

WIpAA ^ ^!(|A| + i^l^N-n-p. 

Also, note that for tt e Lip we have |5®(7r)| = n + p. Then 


n < aa = n i^+<-i n ias..j^(iAi+ivc) 


■-\n+p 


xeS^ (tt) 
TTx =X 


x€S'^(7r) 

■nx=^x 


xeS^{tt) 

We finish the proof by continuing in (13.621) for snfHciently small e and C big enough 


|(-£“/i)^""P^(A) -P‘5(A)| 
^ C{C+\\\Y^'-e°‘ 


+ (|A|+Arc)"^ ^ Ar!(|A| + 

TTGllo i = 0 

N—n N—n—p 

+ YA\+Ncr^^Y E N\{\X\+NcA-^-’’-\e'^hY 

P=1 TTGlIp 1 = 0 


^ + lAD^-^e'" 


GO 

(|A| + Arc)"(A!)2r^(|Al + iVc)'(e“h)' - ll 
^ 1 = 0 

00 

(|A| + mriN - n)(iV!)2[ ^(|A| + iVc)*(£“/r)* - ij 


1=0 

^ C{C + |A|)"”^e“ + 2(|A| + NcY+Ymfe'h + 2(|A| + NcY+Yn - n)(iV!)^£“/i 
^ C(C + |A|)"+^£“. 


(3.66) 


(3.67) 


(3.68) 


□ 
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4 Eigenvectors. Proof of Theorem 12.4 

For the investigation of the eigenvectors we will use some facts presented in m where the author 
considered the spectrum of the unsealed matrix In — P®. Let Xi = —e^hXl, 1 ^ i ^ N, he the 
eigenvalues of this matrix. Of course, the eigenspaces of Q® and Ijv — P® coincide and Theorem [23] (i) 
and Proposition 13.41 imply that 

\\ =0, Ai = 0{e°‘), 2 ^ i ^ n, and ReAf ^ 0, 1 ^ ^ n. (4.1) 

The proof of Theorem 12.41 needs some preparations. For a given subset 7 cz <S consider the first hitting 
time of 7, 


aj := inf{7 ^0: Zl e I}, 
and, for convenience, the first return time to 7, 

t! := inf{7 ^ 1: G 7}. 

For u e C, x e S and I, J c S consider the conditional Laplace transforms 

GT^j{u) := 

and 


(4.2) 


(4.3) 


(4.4) 


(4.5) 


These expressions are finite for every u e C such that Reu ^ 0 but they can be infinite if the real part 
of u is positive. We will discuss the finiteness of these Laplace transforms for certain special cases later. 
Since for any initial state x G S\I we have the equality af = rf, it follows that 


x^IuJ, 

KT;A^) = \i/ xel, 

0, xe J\I. 


(4.6) 


In the following lemma we will prove another useful relation between these two functions (see Section 2 
in [20]). 

Lemma 4.1. Fix an arbitrary x e S, subsets I,JcS and u e C such that Gj’j(u) is finite. Then 


ysS 

Proof. An application of the strong Markov property yields 

GIAA = E,, = J] E,, = y] = y) 

yeS 

yeS 


(4.7) 


(4.8) 


□ 


Now let us introduce the following quantities. Define g N to be the maximal number of steps 
the deterministic motion Z° needs to reach M when starting somewhere in S, 

ymax _ niaxmin{fc 5= 1: Z^x) e M}. (4.9) 

x€S 

Recall the definitions of the interval boundaries Oy and by of ly and of the state y*{x) for a given x e S 
and define 

D := D{h,S) := mindist ^a; — hU'{x), [ay*{a;),j. (4-10) 
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This quantity gives the minimal distance which the deterministic motion has to the interval boundaries 
after one time step and is crucial when it comes to comparing the trajectories of Z° and Z®. In particular, 
(13.51) implies that 

minP4Z? ^ Z!) = 1 - P4Z? = Zl) = 1 - = y*{x)) = ^ (4.11) 

T£<S ^ ' 


We denote p® := £°‘hD~°‘, p® = 0{£°‘) as e —> 0. 

We are interested in the radius of convergence of the Laplace transforms G^’^j and Kj'j in dependence 
on e and for the case 7, J c A4. For this let 

u%i sup |u ^ 0: G^j^{u) = < oo for all s e 5|. (4.12) 

Obviously, for u ^ 0 and I Ai one always has G^’^{u) ^ G^j^(u) so we can immediately conclude 
Gj’^(u) < CO for every u e C with 0 ^ Rew < u%^. 

Proposition 4.2. We have 


lim = 00 . (4.13) 

e—>0 

Proof. Consider fc 5= 1 such that i ^ ^ ^ The occurrence of the event = fc} implies 

that the Markov chain Z® did not follow the deterministic motion Z° at at least one step on the time 
interval [0, fe], namely 

{rXi = fc} c {Zf jc Zi for at least one 1 ^ k}. (4.14) 


Similarly, for + 1 ^ fc ^ gymax occurrence of = k} implies that such a deviation happened 

at least twice and so forth. Hence, for p G N and + 1 < fc ^ (p + 1 ) 7 ’“'“^’= one has 

{tXi = fc} c {Z® ^ zf for at least p time instants {Zi,..., Z^} cz {1,..., k}}. (4.15) 


Therefore for all a; G <S, p G N and pT™®'* + 1 ^ fc ^ (p + 
property that 

Pa;(rXi = k) (^Zf jC zf for at least p time instants 
^ Pj; I Z® jC zf for at least p time instants 


^ s 

l = fL 


ip- 


1)T" 

Z 


iP^ 


l)Tm®x .^g with the help of the Markov 

{Zi,.. 

•,^m} c {1 ,.. 

.,fc}) 


{Zi,.. 

•,^m} c {1 ,.. 

.,(p + l)T”^='}) 

(4.16) 


We apply the well known estimate for the binomial coefficient (^) ^ (x)*^’ n ^ 1, 1 ^ Zc ^ n, and the 
bound p ^ I ^ {p + 1)7’™^’' to get the estimate 


ip- 


1)T" 

Z 




(4.17) 


for some constants Ci > 0 and 77 > 1. Then for some constant C 2 > 0: 

= Zc) ^ C'i((m + l)r““' -p+ l)77''(p®)'‘ ^ C2PiKp^Y. (4.18) 

_^mmax , 

For any ^ 0 we set Cz = X!fc=i = k) and obtain 

00 

^ e^'P^rX, = k) 

fc -1 

OO 

= ^3+S S e“''P4rX, = fc) (4.19) 

M=1 

00 

^ Cs + 774 y; 

This series converges if e“^ 77p® < 1 or, equivalently, if w < yXax In(^) “ ln(77). The claim of the 

proposition follows from the fact that p® is of order 0{£°‘). □ 
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From this proposition it follows that for a fixed uo > 0 one can always choose e small enough such 
that Uo ^ u%^ and therefore we can immediately deduce the following corollary that will provide us with 
an approximation of G^’j{u) with respect to u in a neighbourhood of 0 . 


Corollary 4.3. Fix uq > 0 and eo > 0 such that uq ^ u%i and 0 < e < eq. Let x e S and I c _M, Then 
there exists a constant L > 0, which can be chosen independently of e and x, such that for every u ^ uq 




(4.20) 


Proof. The choice of e guarantees that G^’‘^{u) is finite for u ^ uq. We then use the Taylor expansion 
and obtain 




l| d”G^:M( 0 ) | 
n! I dw" I 


(4.21) 


We have to estimate the expressions - 


| d”G^:X,( 0 ) | 

I du" I 


^ , An 


E |_d_ uk 

(^lin 


k = l 
00 


in the limit of small e. An easy calculation yields 
• Pa:(T/ = fc ^ rXl) 


(4.22) 


k=l 
00 

fe-l 

00 00 

= E 

k=\ l=k 

We proceed with similar arguments and the same constants G 2 and K as in the proof of Proposition 14.21 
For a given fc 5 = 1 let G N be defined by the relation + 1 A: ^ (/r^ + Then, 

aO 00 00 GO (/4+1)T*^^^ 

y r y P4rXi = 0 s; E E E = o 

fc=l Z = ^T™^^4-1 

00 CO 

E fc" E C2KKp^r 

k=l 

00 00 

^ Cs E E 

k=l 
C 3 


k=l l=k 


(4.23) 


1 - 2Kp‘ 


E k"{2Kp^ 


k = l 


We note that fik ^ , abbreviate r® := (2Jfp®)T™^, choose e small enough such that ree“° < 1/2 

and find yet another constant G 4 > 0 such that 

00 I m 00 00 

|G/;Xt(u) - Gium ^ G4 E ^ E E (r^)''(e'='“' - 1) 


= C4 


7l\ 

n=l ‘ k=0 


(el“l - 1 ) 


(4.24) 


(1 - r®)(l - r®el“l) 


^ L|u| 


with a certain constant L > 0. 


□ 


Lemma 4.4. Let A be a nxn matrix such that 0 G (t(A) and eigenspace corresponding to the eigenvalue 
0 be one-dimensional. Let tp be an eigenvector normalized such that ||r/>||oo := maxi 5 jjs;„ = 1. 
Furthermore, let p{e), e > 0, be a sequence of vectors such that lime^o ||/o(£)||oo = 0 and assume that 
for every e > 0 the equation Ax = p{e) possesses a solution x = '0 (e) such that ||0(£)||oo = 1. Then 
lime^o 110 (e) - 0 IIOO = 0. 
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Proof. The matrix A defines a linear mapping from M" to R". Let -Ei := ker(A) and E 2 im(A) 
such that R" = Ei 0 E 2 . Then it is easy to see that the restriction A 2 := A|_b 2 of A on E 2 a bijection 
from E 2 onto E 2 . Indeed, since every x G R'^ such that Ax = 0 is an element of Ei we can deduce that 
ker(A 2 ) = 0 and hence A is injective. We also know that for every x G E 2 , we can hnd y = yi + 1 / 2 , 
yi G El, y 2 e E 2 such that x = Ay = A 1/2 = A 2 y 2 which implies the surjectivity. 

Now let '0(e) be the normalized solution to A0(e) = p{e)- We then can write 0(c) = t{e) + u{s) 
for certain t(e) G Ei and u(c) G E 2 . Hence Au(e) = p(e). Since A is invertible on E 2 and (A) ^|b2 = 
(A 2 )[b 2 we can write u(e) = (A 2 )”^p(c) and therefore ||0(c) — t(c)||oo = ||u(e)||oo —> 0 as e ^ 0. Since 
t(e) e El = {aif : a g R} we can conclude by the normalization condition on 0(e) that ||0(e) — 0||oo ^ 0 
what finishes the proof. □ 

Now we are able to finish the proof of Theorem m 
Proof of Theorem \2.4\ 

For i 2 n, let nf be defined by Af = 1 — e““*. By SHI it follows that lime^o uf = 0. We use a 
convenient representation of the eigenvectors of Ijv — P® obtained in [M Lemma 4.1]. Namely, for any 
c > 0, applying this lemma with I = 0, J = M we get that there is a non-zero vector (0mi, ■ • ■, 0m^) 
such that the eigenvector of liv — P® is expressed in terms of the Laplace transforms as 

n 

0S’' = E («f), * G 5. (4.25) 

j=i 

For definiteness let us normalize maxi 5 jjs;„ |0m*| = L Expanding the functions u K{u) with the help 
of Corollary 14 .3 1 we get 


'1, 

X ^ 

ruj. 




- IPx 



+ bx,mj (e), 

X G flj, 

X A rrij, 




+ 53:,mfc(e), 

X G 7 

k A j, 


where maxa,,*, | 6 a;,mfc(e)| = C>(e“). For x G flj we also have 

liiri Pj; ^(T®^ =1 as well as lirnPaj^a®^ ^ = 0 , fc ^ j. 


(4.26) 


(4.27) 


Recall that 0®’“ satishes the eigenvalue equation A®(Af)0®’* = 0. From this system, let us single out n 
equations indexed by mi,..., mn: 

n 

^ 0®’; 2 (A®(A®))„„, • (P, (a® ^ ^ (e)) = 0 , fc = 1 ,..., n. (4.28) 

i=l ysS 

Abbreviate A'^(A) := Q — Ain and recall the following limits: 

= (A‘^(A))i,j, j = l,...,n, 

(4 29) 

lim e“/iAf = X^. 

Thus we can denote 

n 

Pk{£) = Yj fc = 1,... ,n. (4.30) 

1=1 

where max* |p),(c)j —* 0 as e —> 0. Recall that the eigenspace of X^ is one dimensional and apply Lemma 
m That proves the assertion of the theorem. □ 
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